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ABSTRACT 

Resolved submillimeter imaging of transitional disks is increasingly revealing the complexity of disk 
structure. Here we present the first high-resolution submillimeter image of a recently identified tran¬ 
sitional disk around IRAS 04125+2902 in the Taurus star-forming region. We measure an inner disk 
hole of ^20 AU around IRAS 04125+2902 by simultaneously modeling new 880 /rm Submillimeter 
Array (SMA) data along with an existing spectral energy distribution supplemented by new Discov¬ 
ery Channel Telescope (DCT) photometry. We also constrain the outer radius of the dust disk in 
IRAS 04125+2902 to ~50-60 AU. Such a small dust disk could be attributed to initial formation 
conditions, outward truncation by an unseen companion, or dust evolution in the disk. Notably, the 
dust distribution of IRAS 04125+2902 resembles a narrow ring (Ai?^35 AU) composed of large dust 
grains at the location of the disk wall. Such narrow dust rings are also seen in other transitional disks 
and may be evidence of dust trapping in pressure bumps, possibly produced by planetary compan¬ 
ions. More sensitive submillimeter observations of the gas are necessary to further probe the physical 
mechanisms at work in shaping the spatial distribution of large dust in this disk. Interestingly, the 
IRAS 04125+2902 disk is significantly fainter than other transitional disks that have been resolved 
at submillimeter wavelengths, hinting that more objects with large disk holes may exist at the faint 
end of the submillimeter luminosity distribution that await detection with more sensitive imaging 
telescopes. 

Subject headings: planets and satellites: formation — protoplanetary disks — planet-disk interac¬ 
tions — stars: pre-main sequence — infrared: planetary systems — submillimeter: 
planetary systems 


1. INTRODUCTION 

Transitional disks have been singled out for their pos¬ 
sible connection to planet formation. These disks have 
significant depletions of dust in their inner regions and 
many researchers have proposed that suc h clearing is due 
to pla net-disk interactions (see review bv lEspaillat et all 
1201411 . Potential pro toplanets have also been identified 
in these disks (e.g., lKraus et al.l 120111: Huelamo et al.l 
I2011I iReggiani et al.l D0l4 iBiller et al.l 2014( 1. lending 
further support to the idea that planets are leading to 
the observed disk clearings. How ever, alternative mecha¬ 
nisms such as phot oevaporation (iHollenbach et all 1 1994 
IClarke e t al.l 1200111 can explain a subset of transitiona l 
disks ( Alexander fc Armitagel l200^ lOwen et al.l [201111 . 
A more detailed look at these objects is warranted to 
understand the links to planet formation. 

The detail in which we can study transitional disk 
structure has greatly improved over time. The holes in 
transitional disks were first inferred based on “dips” seen 
in the ground- based infrared (IRl spectral energy distr i- 
butions lSEDs: [Strom et aTIllbsot fSkrutskie et al.llTo^ . 
Spitzer later greatly e xpanded the known number of 
transitional disks ('e.g.. iMuzerolle et al.ll2010l: iKim et al.l 
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1201311 and provided much more detailed SEDs. Spitzer 
IRS a llowed for detailed modeling (e.g., iD’Alessio et al.l 
[200^ . as well as the detection of gaps in disks 
(as opposed to holes ) later dubbed “pre-transitional” 
(|Espaillat et al.l [200711 . Such modeling inferred that the 
gaps and holes in these pre-transitional disks and transi- 
tional disks were very large, of the ord e r of 10s of AU 
(iCaNet et al.l 120051 : iBrown et al.l l2007t lEsoaillat et al.l 
1200711 . Submillimeter images confirmed these predic¬ 
tions by revealing large cavities in the distribution 
of large grains in the disk (e.g.. iHughes et al.l 2007 


Brown etal.ll2(I08l:lHughes et al.ll2009HBrown et ai.l2009 : 


ig. _ 

Andrews et al.l l2009t ~ ll^lla et al.l 1201011 : [Andrews et al.l 

2011hll . 


Resolved submillimeter interferometric imaging of 
disks is greatly expanding our understanding of the spa¬ 
tial distributions of dust and gas. Submillimeter images 
have revealed that the large dust grains in the disk do 


Isella et al.[[2007l: [Panic et al.[[2009l: [Andrews et al.j 

2012 

Rosenfeld et al.l[2013[: de Gregorio-Monsalvo et al.l 

2013 

Pineda et al.l[20l4 [Pietu et alJ[2014[l. The small dust 


disks seen in these observations have been interpreted 
as evidenc e of dust evolution, in particular dust radial 
drift f e.g.. iWeidenschillinel 119771: iBarriere-Fouchet et al.l 


120051 : iLaibe et al.l l2008fb iBirnstiel fc AndrewsI ( 201411 
showed that dust radial drift will lead to the concen¬ 
tration of larger grains at smaller radii, as well as a 
sharp outward trunc ation of the disk, as is observed (e.g., 
[Andrews et al.l[2012ll . This differentiates the signatures 
of dust radial drift from those of dust grain growth, 
which should also lead to the concentration of larger 
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grains towards the inner radii the disk, but not a sharp 
outw ard truncation of the disk (iBirnstiel fc AndrewsI 
[Ml . In the case of transitional disks, dust radial 
drift can lead to concentrations of dust at the outer 
wall as the dust gets trapped due to pressure bumps 
in the gas produced by planetary companions (e.g., 
iZhu et ahl 120121 : iPinilla et al.l l20ll and recent obser¬ 
vations point t o the detection of su c h “dust traps” 
in disks l e-g.. iRosenfeld et al.l 120131: iCasassus et all 


20131: iFukagawa et al. 2(iiT^ van der Marel et alJ 201.1l: 


Isella et al .112 01 311 Perez etMl20l4 iPineda et al.ll20ll !^ 

Note that most of the above studies have focused on 
the brightest submillimeter disks. Submillimeter obser¬ 
vations to date have found that there is a high frequency 
of large cavities amongst the more lumi nous half of the 
disk population. [Andrews et al.l (|2011b[ l estimated that 
at least 1 in 5 of the disks in the bright half of the 
Taurus and Ophiuchus submillimeter luminosity distri¬ 
bution have cavities with radii >15 AU. This suggests 
holes and gaps tend to be found in more massive disks, 
which was also found in a demographic study of Tau¬ 
rus by INaiita et al.l (|2007f) . However, this may be due 
to a selection effect given that few faint disks have been 
resolved. 

Here we add observations of a faint (<20 mJy at 
1.3 mm; [Andrews et al.l[2013f) transitional disk in Tau¬ 
rus to the growing literature on (pre-)transitional disks. 
We present an SMA high-resolution 880 ^m image of 
IRAS 04125-1-2902 as well as new UVBRI photometry 
from the DCT. IRAS 04125-1-29 02 was first identified as a 
candidate member of Taurus bv iKenvon et al.l ([199411 and 
later confirmed as a member bv iLuhman et alT i 20091 ) 
who identified it as an M-type star. A transitional disk 
was inferred usinR Spitz er IRAC and MIPS photome- 
trv (iLuhman et ah 12009( 1 and later Spitzer IRS spectra 
(jFurlan et al.ll201~ ). IRAS 04125-1-2902 joins the pre¬ 
viously imaged Taurus (pre-)transitional disks around 
GM Aur, DM Tan, UX Tau A LkCa 15, and RY Tan 
(iPietu et al.1120061: [Hughes et al.1120091: llsella et al.l[2010a[ : 


[Andrews et al.ll2011blla[: llsella et al.ll2012ll2014fl. 

This paper is organized into the following sections. 
In Section 2, we present the SMA and DCT observa¬ 
tions of IRAS 04125-1-2902. In Section 3, we constrain 
the stellar properties and derive an accretion rate for 
IRAS 04125-1-2902. We also report the results of simulta¬ 
neous broad-band SED and submillimeter visibility mod¬ 
eling. In Section 4, we discuss our results in light of pre¬ 
vious transitional disk studies. Finally, in Section 5, we 
summarize the key results of this work. 


2. OBSERVATIONS AND DATA REDUCTION 
2.1. Submillimeter Data 

IRAS 04125-1-2902 was observed with the 8-element 
SMA interferometer on Mauna Kea, Hawaii using the 345 
GHz dual-sideband receivers in the compact (~10-70 m 
baselines; 2011 Oct 26), extended (~30-230 m baselines; 
2011 Aug 19), and very extended (~70-500 m baselines; 
2011 Sep 8) configurations. The SMA correlator was 
configured to process 4 GHz of bandwidth per sideband, 
centered ±6 GHz from the local oscillator frequency of 
340.8 GHz (880 p-ui), divided into 48 spectral chunks 
containing 32 coarse resolution (3.4 MHz) channels each. 

One chunk centered on the GO J=3—2 transition at 


345.796 GHz had slightly higher resolution (correspond¬ 
ing to a velocity resolution of 0.7 km s“^). Note that 
we did not detect CO emission from our target. We 
do not derive an upper limit since the line is not de¬ 
tected and there is substantial molecular cloud contam¬ 
ination in the CO line. The observations cycled be¬ 
tween IRAS 04125-1-2902 and the nearby quasars 3C 111 
(1.8 Jy) and J0510-I-180 (0.9 Jy) over ^10-15 minute 
intervals. Bright quasars (3C 84, 3C 279) and primary 
flux calibration sources (Uranus, Callisto) were observed 
when the target had low elevation. Observing conditions 
were generally good, with precipitable water vapor levels 
under 2 mm. The raw visibilities were calibrated us¬ 
ing the MIF0 software package. The passband response 
was corrected based on observations of bright quasars, 
and the amplitude scale was referenced to Uranus and 
Callisto. Complex gain calibration was performed using 
the repeated observations of the nearby quasars. The 
visibilities were spectrally binned into wideband contin¬ 
uum channels; the data from all three configurations were 
then combined. The calibrated visibilities were Fourier 
inverted (using a Briggs weighting with robust param¬ 
eter of -1, chosen by experimentation to represent the 
best combination of angular resolution and S/N for these 
data.), deconvolved, and then restored with a synthe¬ 
sized CLEAN beam using standard tasks in the MIRIAD 
package. The resulting image, shown in Figure 1, has a 
synthesized beam size of 0733 x 0726, and an RMS noise 
level of 1.3 mJy beam“^. We measure a flux density 
of 0.027±0.001 Jy for IRAS 04125-L2902. The flux den¬ 
sity was estimated from the shortest visibility spacings, 
but one finds the same values (within the uncertainties) 
from integrating the synthesized images within the Icr 
contours or a reasonable aperture (170 in radius). 

2.2. Optical Data 

On 2013 Nov 29 we used the Large Monolithic Imagei0 
(LMI) on the 4.3 m DCT to obtain UBVRI photometry 
of IRAS 04125-L2902. The FOV of the LMI is 12f5 x 
12f5 (0712 per unbinned pixel). We used 2x2 pixel bin¬ 
ning for these observations (0724 per pixel). The bias, 
flat-held, and overscan calibration of the CCD images 
were performed using the ccdred package in IRAF. The 
photometric calibration of all images was carried out in¬ 
teractively using the daophot and phot cal packages in 
IRAF, with standard stars selected fro m [Landoltl (Il992() 
and fo llowing; the procedure outlined in iMassev &: Da\^ 
()1992l l. Photometry is presented in Table 1 and Figure 2. 
Note that our data of IRAS 04125-1-2902 are consistent 
with previous B- (17.83, 17.44), R- (14.39, 14.34), and 
/-band (12.81) measurements from the US Naval Obser¬ 
vatory catalog (|Zacharias et al.l [2009D given the obser¬ 
vational and calibration uncertainties. For clarity, we 
do not show these additional data in Figure 2. Our 
IRAS 04125-1-2902 data are also consistent with previ¬ 
ous u grz SDSS data (Figure 2: IAdelman-McGarthv et ahl 
[2011(1 within the uncertainties. 

3. ANALYSIS & RESULTS 

In this section, we first review the disk model used 
in this work. We then discuss the adopted and derived 

® http://www.cfa.harvard.edu/ cqi/mircook.html 

® http://www.lowell.edu/techSpecs/LMI/LMI.html 
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Figure 1. SMA 880 fim continuum emission of IRAS 04125+2902. From the left, panels show the SMA continuum image, the synthetic 
model image, the imaged residual, and the visibilities (black points) with the model visibility (red line). Counter intervals are 3 cr. There 
is a double-peaked emission structure and a null in the visibilities, indicating a large cavity in the disk. 


Table 1 

Optical Photometry 


'i'arget 

U 


V 

R 

I 

IRAS 04125-1-2902 

19.191 ± 0.041 

17.442 ± 0.011 

15.487 ± 0.003 

13.954 ± 0.005 

12.418 ± 0.005 


stellar properties and accretion rate, which are used as 
inputs in the model. Finally, we present the first detailed 
modeling of IRAS 04125+2902. In short, we generated 
synthetic images for a grid of disk models with the best 
chi-squared values from SED fitting and compared them 
to the observed submm images to arrive at our best fit. 
Individually, SED modeling and submm image modeling 
have their limitations, but below we discuss how mod¬ 
eling SEDs and submm images together can help break 
some model degeneracies. 

3.1. Disk Model 

He r e we use the disk models of ID’Alessio et al.l (I1998L 
11999112001LI2005LI2006D . The application of these mod¬ 
els to the particular case of fittin g transitional disks ha s 
been discussed in detail before bv lEspaillat et al.l (I2010D . 
In short, the D’Alessio et al. models are of irradiated, ac¬ 
cretion disks around pre-main-sequence stars. The tem¬ 
perature and density structure of the disk is calculated it¬ 
eratively and the surface density of the disk is tied to the 
accretion rate. A transitio nal disk model has a hole rel¬ 
ative to a full disk model (|D’Alessio et al.l[2006D . There 
is a frontally illuminated “wall” at the edge of the hole 
(i.e., at the inner edge of the disk). This wall dominates 
the mid-infrared (MIR) emission (from ^20-30 ^m) and 
the disk dominates the emission beyond '^40 /rm. 

The composition of dust used in the disk model im- 
pact s the resulting emiss ion and derived disk properties 
(see lEspaillat et al1l201ClD . Here we include silicates and 
graphite with fractional a bundances of 0.004 a nd 0.0025, 
respectively, following the lDraine fc Led (II984I) model for 
the diffuse ISM. We calculate the silicate (olivine) and 
grap hite opacities using Mie theory and optical constant s 
from iDorschner et al.l (lIOO'sl) and iDraine fc Led (|I984D . 
respectively. The models assume spherical grains with a 
size distribution that scales like a~^ between grain radii 
of amin and Umax and p is 3.5 (|Mathis et al.lTl97^ . In 
the wall and the disk atmosphere, amin is held fixed at 
0.005 /im while Omax is varied to achieve the best fit 
to the SED, particularly the 10 pm and 20 pm silicate 
emission features. Near the disk midplane, the maximu m 
grain size is held fixed at 1 mm (|D’Alessio et al.l[2006[) . 

In fitting the outer wall of the disk, Twaii, the tempera¬ 


ture at the surface of the optically thin wall atmosphere, 
is adjusted to fit the SED. The height of the wall {zwaii) 
is equal to the disk scale height (H) at the the radius in 
the disk at which the wall is located (Rwaii)- Rwaii is 
derived using the best fitting Tyjaii following 


i?. 


wall 


-|- Lace) ^2 -|- \ 

lOTTCTfi Kd 


-I 1/2 


T/ 


( 1 ) 


'all 


where an is the Stefan-Boltzmann constant, Td is the 
total mean optical depth in the disk, and Kg and Kd 
are the mean opacities to the incident and local radi¬ 
ation, respectively. L* is the stellar luminosity and Lace 
(^ GM^M/R*) is the luminosity of the stellar accretion 
shock. 

The parameters of the disk which are varied to fit the 
data are the viscosity parameter (a) and the settling pa¬ 
rameter (e). One can interpret varying a for disks with 
the same inner and outer radii as fitting for the disk mass 
since M^isfe oc M/a (see Equation 38 in iD’Alessio et al.l 
II998D . We note that the mass accretion rate onto the star 
does not necessarily reflect the mass transport across the 
disk nor do we expect that the mass accretion rate is con¬ 
stant throughout the disk. However, for simplicity, here 
we assume that the mass accretion rate measured onto 
the star is representative of the disk’s accretion rate. We 
parameterize settling as a depletion of dust in the upper 
layers using e = Catm/Cstd (i-e., the dust-to-gas mass ratio 
of the disk atmosphere divided by the standard value), 
with a corresponding increment of the dust-to-gas mass 
ratio near the midplane. One can also vary the disk 
inclination and outer radius (Rwaii), which can be con- 
strained with high-res olution submillimeter images (e.g., 
[Andrews et ahll^Ilbl) . Details of our fitting procedure 
in the particular case of IRAS 04125+2902 can be found 
in Section 


3.2. Stellar Properties 

The star is a signiheant source of heating for the disk 
and so the stellar properties are important input param¬ 
eters for the code. The stellar parameters used here 
are listed in Table 2. We use the same spectr a l type 
for IRAS 04125+2902 adopted in iFurlan et al.l (|2Cill[) 
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Figure 2. Disk model (black solid line) fit to the dereddened 
SED (blue) of IRAS 04125+2902. The original data (red) were 
dereddened (blue) with the extinction listed in Table 1 and de¬ 
tails can be found in Section 3.2. Refer to Section 3.4 and Table 2 
for model details. Here we show new DCT data (open circles; 
Table 1) from this work and S MA (open squares) data first pre¬ 
sented Wealso include SDSS data (open 

stars: |A delman-McCar thv et al.ll2QlUh JHK data (closed circles ; 
ISkrutsk'i^^^aLil2006l) . WISE (closed triangles; [Wrigh^^ay[2^R|), 
IRAC and MIPS p hotometry (ope n triangles; ILuhma n et al.ll2010l : 
IRebull et al.lf^iol'l. IRAS ('crossesIBeichman e^ahl^QSSlh AKARI 

(close d sauar ^Ishihara et al.l2010ll , and IRS spectra llFurlan et al.l 

I201U') . The broken ma genta line represents the stel lar photosphere 
based on colors from IKenvon HartmannI l| 199511 and scaled at 
J. Separate model components are the stellar photosphere (ma¬ 
genta dotted line), the outer wall (red dot-short-dashed line), and 
the disk (cyan dot-dashed line). (A color version of this figure is 
available in the online journal.) 


and [And rews et 'HI ^201,Sfl . which originally conies from 
iLuhman et al.l ( 2009tl using IR spectra and with reported 
uncertainties of ±0.25 subclasses. We measured a vi¬ 
sual extinction by c omparing observed co l ors to pho- 
tospheric colors from IKenvon fc HartmannI (119951 ). We 
adopted an Ry of 3.1 (Ay/Aj=3.55) and dereddened 
our data using the extinction curve from lMcClurel (l2009f ) 
for IRAS 04125±2902 {Av< 3). Note that our derived 
extinction value for IRAS 0 4125±2902 (A \/=2.7±0.5) 
agrees with that fou nd by iFurlan et al.l (I2011D and 
[Andrews et al l (120131). The stella r temp erature (T*) 
is taken from (Kenyon fc HartmannI (119951) . based upon 
the adopted spectral type. The luminosity is calcu- 
lated with dereddened J-band pho tometry from 2MASS 
( Skrutskie et al.l l2006f ) following IKenvon fc HartmannI 
( 1995D assuming a distance of 140±20 pc for Taurus 
( Bertout et al.lll9^ . R* is calculated using the derived 
luminosi ty and adopted t emperature. M* was derived 
with the lSiess et al.l (|2000f ) pre-main sequence evolution¬ 
ary tracks using T* and L*. The uncertainties in L* and 
R* are dominated by the uncertainty of the extinction 
measurement. We note that our values for the luminos¬ 


ity and_jnuss_ajeingood agreement with those found 
bv [Andrews et al.l (I2013D with a Bayesian inference ap¬ 
proach. 


3.3. Accretion Rate 

In the current paradigm of TTS accretion, material 
from the disk is channeled onto the stellar surface via the 


Table 2 

Stellar and Disk Properties of 
IRAS 04125+2902 


Stellar Properties 

M* (Mq)^ 

O.SiO.l 

R* (R©)! 

l.SiO.l 

T* (K)2 

3720±70 

L. (L©)! 

0.40±0.06 

Av' 

2.7±0.5 

Spectral Type^ 

M1.25±0.25 

Inclination (°)® 

<30 

M (M© yr-l)l 

31® X10-1° 

-Outer Wall ± Disk- 

^max (MI^)'^ 

10-100 

T'n,aU (K)3 

85-95 

(AU)4 

1.8-2.5 

R»aH (AU)4 

18-24 

R'disfc (AU)^ 

50-60 


0.01-1.0 


0.0001-0.0003 

Mdisfc (M©)'* 

0.005-0.007 


^ These values were derived in this 
work; see Sections 3.2 and 3.3 for 
more details. 


^ The spectral type was adopted 
from ILuhman et al.l 1 120091) and 
the stellar temperature is from 
IKenvon Sz Hartman 3 using 

this adopted spectral type. 

^ These values are free parameters 
in the model and were derived based 
on fitting the SED and submm im¬ 
age; see Section 3.4 for more details. 

The wall height, wall radius, and 
disk mass were calculated based on 
the best-fitting disk parameters fol¬ 
lowing the methods discussed in 
Section 3.1. 


stelly magnetospheric field lines le.g.. lUchida fc Shibatal 
I984I). forming a hot acc retion shock (^5,000-10,000 K 
Calvet fc Gullbriii^ 1199811 . Accretion rates can be de¬ 
rived by measuring the accretion luminosity due to the 
accretion shock on the stellar surface, which emits excess 
emission ab ove the stellar photospher e, that peaks in the 
ultraviolet (iCalvet fc Gullbrin3ll998D . 17-band excesses 
have been fo und to correlate with the total accretion 
shock excess (|Calvet fc Gullbriii^ 1199811 and so [/-band 
data can measure the ultraviolet excess above the stellar 


photosphere (e.g.. iGullbrine et alJ 119981: 

White & Ched 

120011 IRigliaco et al.1120121: IFsoaillat et al. 

I 20 TT 


Here we measure an accretion rate for 


IRAS 04125T2902 of 3xl0-i° Mq yr"! using the 
DCT [/-ba nd data presen t ed in this work and the 
relation in iGullbring et al.l (I1998D . This method has 
a typ ical uncertainty of a factor of 3 (| Calvet et al.l 
I2004f) . This low accretion rate is consistent with 
IRAS 04125±29 02’s previously repor ted Ha equivalent 
width of 2.3±0.3 (|Luhman et al.ll2009D. which is expecte d 
from a non-accreting, weak-TTS (|White &: Basrill2003D . 
Note that at such low M, [/-band emission may be 
domi n ated by chromospheric emission (jHoiidebine et aP 
1199(11 : iFranchini et al.l 119981) , which can contaminat e 
mass accretion rate measurements (|Inglebv et al.l[20Tll) . 
Future high-resolution ultraviolet spectra are needed 
to reliably constrain the accretion rate of this object 
further. 
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3.4. BED and Submillimeter Image Modeling 

Using the stellar properties in Table 2, we ran a grid 
of '^2,300 models where we varied s-max (5 /im, 10 ^m, 
100 ^m), T^aii (75 K, 85 K, 95 K, 105 K), e (0.01, 0.1, 
0.5, 1.0), a (0.01, 0.001, 0.0001, 0.0003), Rdisfc (40 AU, 
50 AU, 60 AU, 100 AU), and the inclination (20°, 40°, 
60°). The parameters explored for amax, Tiwaii, and 
R-disfe were constrained based on inspection of the SED 
and submm image. The lack of strong silicate emission 
features in the IRS spectrum points to large (>5 /rm) 
grains in the disk atmosphere. The steep slope of the 
IRS spectrum longwards of 10 ^m suggests cooler tem¬ 
peratures for the wall atmosphere (>100 K). The submm 
image points to a relatively small (<100 AU) disk. The 
parameters explored for e, a, and the inclination were 
chosen to cover a broad range of values with reasonable 
sampling while avoiding creating a computationally pro¬ 
hibitive grid. We note that the above grid is not meant 
to be exhaustive or lead to a unique fit. We aim to iden¬ 
tify the best fits to our data from a reasonable parameter 
space meant to guide future observations. 

For the models with the best chi-squared fits to the 
SED, high resolution model images at 880 /im were gen¬ 
erated, and their Eourier transforms were sampled at the 
same spatial frequencies as the SMA observations. These 
model visibilities were then compared with the data in 
the for m of an azimuthallv-a veraged deprojected visibili¬ 
ties Isee lAndrews et ^1120091 1. as shown in the right-hand 
panel of Eigure 1. Synthesis images were generated from 
the model and residual visibilities in the same way as for 
the data (middle panel of Figures 1). 

A representative best fitting model to the submillime¬ 
ter image and SED is shown in Figures 1-2, respec¬ 
tively. The particular model shown has Ru,aiz=20 AU, 
Rciisfc=50 AU, e = 0.5, a = 0.0003, *=20°, and 
amax=100 /im. We do not show each of our best-fitting 
models since the fits are very similar. Instead, in Table 2 
we list the parameter ranges for the best fitting mod¬ 
els. The best fits to both the SED and submillimeter 
visibilities of IRAS 04125+2902 were achieved with disk 
models with wall radii of 18-24 AU and disk radii of 50- 
60 AU. We can also constrain the inclination of the disk 
of IRAS 04125+2902 to <30°. (We note that we arrived 
at this constraint by running additional models with in¬ 
clinations between 0° and 40° in increments of 5° for the 
best fitting model shown in the figures.) The wall and 
disk radii derived here as well as the inclinations were 
driven by the fit to the submillimeter visibilities since 
submm/mm imaging is sensitive to the spatial distribu¬ 
tion of large dust in the disk midplane while the SED is 
not as sensitive to the spatial distribution of dust. The 
wall radius is determined from the null position in the 
visibilities. All showed that the submm/mm visibilities 
are sensitive to changes in the wall radius of about 10% 
given typical noise levels. See All for a more in depth 
description regarding the implications of disk ring mor¬ 
phologies and nulls in submm visibilities. 

The derived maximum grain size in the wall and disk 
atmosphere is largely based on fitting the IR SED since 
it is more sensitive to smaller dust grains in the upper 
layers of the disk. In general, the dust grain size in the 
disk atmosphere can be constrained by the silicate emis¬ 
sion features in the Spitzer IRS spectrum. Here we can 


get reasonable fits to the SED with amax between 10 to 
100 /tm. The a and e parameters were co nstrained by fit¬ 
ting t he IR and submm/mm SED data. lEsnaillat et al.l 
(l2012fl found that disk models with the same e-to-a ra¬ 
tio will produce very similar emission in the IR but sub¬ 
stantially different emission in the submm/mm. Hence, 
submm/mm data are necessary in conjunction with IR 
data to break the degeneracy between dust settling and 
disk mass. Since in this work we have both IR and 
submm/mm data, we can break this degeneracy. Here 
we find models with e between 0.01 to 1 and a be¬ 
tween 0.0001 and 0.0003 reproduce the SED and submm 
image best. Note that none of our models reproduce 
the IRAS/AKARI data and we exclude this data from 
our chi-squared fitting. It is unclear with the avail¬ 
able data if the relatively high the IRAS/AKARI fluxes 
can be attributed to variability, as has been seen in 
other (pre-)transitional disks fe.g..lMuzerolle et al.ll2010l 

lEsnaillat et al.l[20Tlt iFlahertv et al.ll2012ll . We leave fur¬ 

ther exploration of this to future work. 

We conclude that there is a relatively narrow 
AR'^35 AU ring of large dust grains around 
IRAS 04125+2902. Our submillimeter images are 
not of high enough sensitivity to confirm or exclude disk 
asymmetries or to test the extent of its dust disk relative 
to the gas disk. We leave this and more detailed studies 
of the dust properties to future multi-wavelength work, 
which would be ideally suited to probe this. 

4. DISCUSSION 

IRAS 04125+2902 is one of the faintest transitional 
disks imaged to date. lAndrews et al.l (l2011bl l found that 
at least 1 in 5 submillimeter-bright, and therefore mas¬ 
sive, disks are (pre-)transitional disks. However, here we 
find a hole in a relatively low submillimeter luminosity 
disk. Given that the lower luminosity end of the disk sub¬ 
millimeter brightness distribution has yet to be probed, 
this suggests that many such disks may contain holes and 
gaps. 

The dust in the disk around IRAS 04125+2902 only 
extends from about 20 AU to 50-60 AU. This is remi¬ 
niscent of a ring-like large dust grain distribution. We 
stress that our SMA observations only trace the submil¬ 
limeter dust continuum and so this derived disk radius 
applies only to the large dust grains in the disk and not 
necessarily the disk overall. Due to the faintness of this 
object, we do not have high quality submillimeter data 
of the gas in the disk and cannot distinguish between a 
disk where both the gas and dust extend over small spa¬ 
tial scales or a disk where the dust is truncated to sig¬ 
nificantly smaller radii than the gas. With this in mind, 
below we speculate on the physical mechanisms that may 
underlie the observed spatial distribution of dust in the 
disk of IRAS 04125+2902. Higher sensitivity observa¬ 
tions of the gas, preferably targeting dense gas tracers 
(e.g., CO isotopologues) given cloud contamination, are 
needed to explore the below scenarios further. 

4.1. Small disk radii 

Future submillimeter observations of the gas in the disk 
of IRAS 04125+2902 may reveal that the gas traces the 
dust distribution. In this case, such small gas and dust 
disks could be attributed to a wide rang e of initial con- 
ditions at the time of disk formation. lAndrews et al.l 
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(|2010f) found an correlation for disks in 

Ophiuchus in which fainter disks are smaller and less 
massive. Even though this work uses a different disk 
model to derive disk masses and outer radii, our re¬ 
sults a re roughly consistent with those of lAndrews et ahl 
(1201011 . Those authors propose that this mass-radius cor¬ 
relation could be the result of a wide range of initial con¬ 
ditions and viscous properties at time of disk formation. 

Alternatively, smaller disks could be the result 
of outward truncation by an unseen companion. 
IRAS 04125-1-2902 is a 4"0 binary, which corresponds 
to a separation of ^560 AU at 140 pc (A. Kraus, 2014, 
private communication). This companion is unlikely to 
be responsible for truncati ng the disk down to 50-60 AU 
according to the models of lArtvmowicz fc Lubowl (|1994[1 
which show that the disk should be truncated to 0.2- 
0.5 a, where a is the binary semi-major axis. However, a 
caveat is that we do not know the eccentricity of the com¬ 
panion’s orbit and disk truncation due to a highly eccen¬ 
tric companion (e >0.75) is not well constrained. There 
are no equal brightness companions outside of '^10 AU 
and no companions near the substellar boundary out¬ 
side '^30 AU (A. Kraus, 2014, private communication). 
Therefore, if the disk of IRAS 04125-1-2902 is being out¬ 
wardly truncated by a companion, this companion would 
have to be <50 Mj^pUer- 

If the gas around IRAS 04125-1-2902 is found to extend 
to much larger radii than the dust in the disk or multi¬ 
wavelength resolved data show a size-wavelength anti¬ 
correlation, then the small dust disk observed here may 
be evidence of dust evolution. Resolved submillimeter 
observations are increasingly revealing that large dust 
grains in disks are distributed on more compact spatial 
scales than the smaller grains and gas in the disk. Eor 
example, the mm-sized dust in the transitional disk of 
TW Hya is truncated at 60 AU while CO gas emission 
is de tected much further out in the disk to about 215 
AU ([Andrews et al.l [201211 . Evidence for compact dust 
emission is also seen in full disks. The full disk around 
AS 209 was observed at four wavelengths between 0.88 
mm and 9.8 mm; at longer wavelengths, the dust disk is 
more compact, a signature that larger grai ns are more 
conce ntrated towards the center of the disk (| Perez et al.l 

IMl). 

Observations revealing concentrations of larger dust 
grains towards the center of the disk have been inter¬ 
preted as evidence of dust evolution. Larger dust grains 
are not tied to the gas and ar e expected to drift radially 
inward (|Weidenschillinsjll977f) . This is due to the differ¬ 
ence in velocities between gas and dust of different sizes. 
Dust travels at the Keplerian velocity while the gas is 
partially pressure supported. In general, the gas pres¬ 
sure decreases with radius so the gas pressure gradient 
is negative and the gas travels at sub-Keplerian veloci¬ 
ties. Smaller dust grains with slower radial velocities are 
coupled to the gas. Larger dust grains experience a head¬ 
wind from the gas; the larger dust slows down and loses 
angular momen tum, spiraling in towa r ds the star. More 
recent work by iBirnstiel fc Andrews! (|2014f l has shown 
that this process will not only lead to more compact dis¬ 
tributions for larger dust grains, but that there should 
also be a distinct, sharp o uter edge in the dus t distri¬ 
bution, as is observed le.g.. [Andrews et 3111201211 . Grain 
growth, in contrast, would lead to a more tapered edge 


(|Birnstiel fc Andrewsll20l3) . 

4.2. Narrow dust rings 

The distribution of the large dust grains in 
IRAS 04125-1-2902 is narrow (from about 20 AU to 50- 
60 AU) and ring-like. As stated earlier, more sensi¬ 
tive observations of the gas are needed to more firmly 
characterize the disk’s spatial properties. However, if 
it is revealed that the dust in IRAS 04125-1-2902 is 
more compact than the gas, IRAS 04125-1-2902 would 
be very similar to other transitional disks where con¬ 
centrations of l arge grains have been detected at the 
outer wall (e.g., iRosenfeld et al.l 1201311 . This has been 
interpreted as evidence of dust trapping due to pres- 
sure bumps in t he gas, possibly i nduced by a planet 
(jZhu et al.l l201^jPjni]la et_aJj 20 In particular, sim¬ 
ulations bv iPinilla et al.l (|2012 ) show that this effect can 
lead to narrow rings of ^20 AU in width at the loca¬ 
tion of the outer wal l . In t he transitional disk of V4046 
Sgr, IRosenfeld et al.l (1201311 measure a ring of 16 AU in 
width at 37 AU. This is five times more compact than 
the gas in this disk, suggesting an accumulation of solids 
at the local gas maximum. In the pre-transitional disk 
of HD 100546, the CO (3-2) gas extends out to 350 AU, 
while a ^ 20 AU ring of larg e dust grains is located at 
^^50 AU ([Pineda et al.l 1201^ . Similarly, narrow rings 
of large dust grains are found in the (pre- ) transitional 
disks SAO 206462. SR 21 (IPerez et al.l[2(^ . and [PZ99] 
J160421.7-213028 ([Zhang et al.ll2014[l . 

The above scenario of a dust trap formed by 
planets is compatible with the interpretation that 
the h oles in disks are being carved out by planets 
fe.g.. iPa ardekooper fc Mellemaj l20Ol iZhu et all 120111 : 
iDodson-Robinson fc Salvkl 1201111 . However, the accre¬ 
tion rate we measure here for IRAS 04125-1-2902 is 
very low (3xl0“^° Mq yr“^), making the hole in this 
disk a possible c andidate for clearing due to photo- 


evaporation (e.g.. iHollenbach et al.l 11994 IClarke et al. 

2001 

; [Alexander et al. 2006 

: Alexander & Armitage 

2007 

; [Alexander et al.ll2013ll. 

Owen et al.l (20111 found 


that they could explain disk holes with radii <20 AU 
and accretion rates <10“® Mq yr“^ with photoevapo- 
rative clearing. However, a low accretion rate is also 
compatible with clearing due to a companion. Be¬ 
sides CoKu Tau/4, IRAS 04125-1-2902 has the lowest 
measured accretion rate of the (pre-)transitional disks 
in Taurus ([Inglebv et al.l l2?)T^ Espaillat et al. 2014). 
The low accretion rate of CoKu Tau/4 (<10“^° Mq 
vr~^: iCohen fc Kuhil 1197911 and hole size (^14 AU; 
iNagel et al.l l20ldir are consistent with clearing due to 
the stellar mass companion detected within the disk 

hole ([Ire land fc Krausll200^IArtvmowicz fc Lubowlll994 

iLubow et al.l 11999(1 . More work is needed to ex¬ 
clude or detect companions within the disk hole of 
IRAS 04125-b2902. 

Based on the above, IRAS 04125-1-2902’s accretion rate 
may indicate its hole is due to disk dispersal by pho¬ 
toevaporation or there is a massive companion in the 
inner disk. The first inte rpretation woul d be c onsis- 
tent with the proposal of lOwen fc Clarkel (1201211 that 
there are two distinct populations of transitional disks 
- those with low mm fluxes, small <10 AU holes, and 
low (<10“® Mq yr“^) accretion rates and those with 
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high mm fluxes, larger holes, and higher accretion rates. 
The former population is consistent w ith models of pho¬ 
toevaporation while th e latter are not (|Owen et al.ll201ll : 
I Owen fc Clarkd 1201211 . The former are also consistent 
with holes in disks developing later in the lifetime of the 
disk. Interestingly, Taurus has both populations. The 
second interpretation of dynamical clearing by a stellar 
or planetary mass companion would need gas observa¬ 
tions to confirm that indeed a narrow ring of dust and gas 
exists at larger radii since, to the best of our knowledge, 
photoevaporation models do not predict such a feature. 


Note that along with narrow rings in (pre-) transitional 
disks, asymmetries have also been detected, which have 
been interpreted as evidence of planet-induced vortices 


(1 Perez et al.l 

2014 ICasassus et al.ll2013HFukaEawa et al.l 

201,31 van c 

er Marel et al.l 1201.31: llsella et al.1 1201.31: 

Pineda et al. 

l2014fl. Our observations are not of high 


enough signal-to-noise to confirm or exclude the presence 
of asymmetries in IRAS 04125+2902. 


5. SUMMARY 

We performed the first detailed characterization of 
the recently identified transitional disk around IRAS 
04125+2902 in Taurus. We presented a new SMA high- 
resolution 880 /rm image and DCT photometry. Using 
simultaneous SED and submillimeter visibilities model¬ 
ing, we find the following: 

1. IRAS 04125+2902 is the faintest transitional disk 
imaged in detail to date, suggesting that there are more 
faint, less massive disks with inner holes that have not 
yet been detected. 

2. IRAS 04125+2902 has a dust disk radius of 50-60 
AU. Gas observations of IRAS 04125+2902 are necessary 
to explore the spatial distribution of the dust relative to 
the gas. If the dust and gas are both truncated to simi¬ 
larly small outer radii, this could point to a wide range 
of initial conditions in disks or an unseen companion. If 
the dust is significantly more compact than the gas in 
the disk, this could be evidence of dust radial drift in 
disks. 

3. IRAS 04125+2902’s dust is distributed as a narrow 
ring (Ai?^35 AU) of large grains at the location of the 
disk wall (^20 AU). Such narrow dust rings are also seen 
in other transitional disks and may be evidence of dust 
trapping in pressure bumps, possibly produced by plane¬ 
tary companions. However, here we also measure a very 
low accretion rate for IRAS 04125+2902 of ^3x10“^° 
Mq yr“^, making this object a potential candidate for 
clearing due to photoevaporation. More sensitive sub¬ 
millimeter observations of the gas in this disk should be 
able to distinguish between clearing due to planets or 
photoevaporation in the disk. 

ALMA will play an important role in probing the lower 
luminosity end of the transitional disk population and re¬ 
vealing the true frequency of small dust disks and narrow 
rings in transitional disks. 
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